Neutral nuclear core vs super charged one 
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Based on the Thomas-Fermi approach, we describe and distinguish the electron 
distributions around extended nuclear cores: (i) in the case that cores are neutral 
for electrons bound by protons inside cores and proton and electron numbers are 
the same; (ii) in the case that super charged cores are bare, electrons (positrons) 
produced by vacuum polarization are bound by (fly into) cores (infinity). 

Equilibrium of electron distribution in neutral cores. In Refs. , the Thomas- Fermi 

approach was used to study the electrostatic equilibrium of electron distributions n e (r) 
around extended nuclear cores, where total proton and electron numbers are the same N p = 
N e . Proton's density n p (r) is constant inside core r < R c and vanishes outside the core 
r > R c , 

n p (r) = n p d(R c - r), (1) 
where R c is the core radius and n p proton density. Degenerate electron density, 

<r) = ^(Pf) 3 , (2) 

where electron Fermi momentum P e F , Fermi-energy S e (Pf) and Coulomb potential energy 
Vcoui(^) are related by, 

£ e (P e F ) = [(P e F c) 2 + my] 1 ' 2 - m e c 2 - V cml (r). (3) 

The electrostatic equilibrium of electron distributions is determined by 

£e(P?) = 0, (4) 

which means the balance of electron's kinetic and potential energies in Eq. ([3]) and degenerate 
electrons occupy energy- levels up to +m e c 2 . Eqs. (12113111) give the relationships: 
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V 2 oul (r)+2m e c 2 V c Ur) ; (5) 
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n e (r) = 37r2( ^ )3 + 2m e cV coul (r)j"' _ . (6) 

The Gauss law leads the following Poisson equation and boundary conditions, 

AV cou i(r) = 47ra [np(r) - n e (r)] ; V CO ui(oo) = 0, V^ui(0) = finite. (7) 

These equations describe a Thomas- Fermi model for neutral nuclear cores, and have numer- 
ically solved together with the empirical formula HQ and /3-equilibrium equation |3( for 
the proton number N p and mass number A — N p + N n , where N n is the neutron number. 

Equilibrium of electron distribution in super charged cores In Ref. HQ 

, assuming that 

super charged cores of proton density §Q are bare, electrons (positrons) produced by vacuum 
polarization fall (fly) into cores (infinity), one studied the equilibrium of electron distribution 
when vacuum polarization process stop. When the proton density is about nuclear density, 
super charged core creates a negative Coulomb potential well —V cou \(r), whose depth is 
much more profound than —m e c 2 (see Fig. [I]), production of electron-positron pairs take 
places, and electrons bound by the core and screen down its charge. Since the phase space 
of negative energy-levels e(p) 

e (p) = [(pc) 2 + m e V] 1 / 2 -Koui(r), (8) 

below —m e c 2 for accommodating electrons is limited, vacuum polarization process com- 
pletely stops when electrons fully occupy all negative energy-levels up to — m e c 2 , even electric 
field is still critical. Therefore an equilibrium of degenerate electron distribution is expected 
when the following condition is satisfied, 

e(p) = [(pc) 2 + my} 1/2 -V coul (r) = -m e c 2 , p = P e F , (9) 

and Fermi-energy 

£ e (P e F ) = e(Pf ) - m e c 2 = -2m e c 2 , (10) 

which is rather different from Eq. (J3J). This equilibrium condition fllUI) leads to electron's 
Fermi-momentum and number- density ([2]), 
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V 2 oul (r) ~ 2m e c 2 V coul (r) ; (11) 



1 



<r) = ^ [Ko» - 2m e c 2 V conl (r) . (12) 
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which have a different sign contracting to Eqs. (I5ll6l) . Eq. (JTj) remains the same. However, 
contracting to the neutrality condition N e = N p and ra e ( r ) Ir^oo — > in the case of neutral 
cores, the total number of electrons is given by 

AP on = / ATir 2 drn e {r) < N p , (13) 
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where tq is the finite radius at which electron distribution n e (r) ( jl"2l vanishes: n e (ro) = 
, i.e., V cou i(r ) = 2m e c 2 , and n e (r) = for the range r > r . N lon < N p indicates that 
such configuration is not neutral. These equations describe a Thomas-Fermi model for 
super charged cores, and have numerically {4] and analytically solved with assumption 
N p = A/2. 

Ultra-relativistic solution In analytical approach 5|, |6) , the ultra-relativistic approximation 
is adopted for V cou \(r) ^> 2m e c 2 , the term 2m e c 2 V COVl \{r) in Eqs. (I5ll6llll|ll2p is neglected. It 
turns out that approximated Thomas-Fermi equations are the same for both cases of neutral 
and charged cores, and solution Vc OU i(?") = hc(37i 2 n P Y^ 3 4>(x) , 
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, for x < 0, 



, . , - • 1 + S i n h(3.44 - V3x) 

nx) = { ^ 1 J } , (14) 

(x+1.89)' x > 



where x = 2(n /3) 1 / 6 'a 1 ■ /2 'n^ 3 (r — R c ) ~ 0.1(r — i? c )/A 7r and the pion Compton length A„- = 
Ujim^c). At the core center r = 0(x — > —00), V CO ui(0) = hc^ir 2 ^) 1 / 3 ~ m n c 2 . On core 
surface r = R c (x = 0), V cou \(R c ) = 3/4V cou i(0) ^> m e c 2 , indicating that the ultra-relativistic 
approximation is applicable for r < R c . This approximation breaks down at r > r . Clearly, 
it is impossible to determine the value r out of ultra-relativistically approximated equation, 
and full Thomas-Fermi equation ([7]) with source terms Eq. ([6]) for the neutral case, and 
Eq. (|T2|) for the charged case have to be solved. 

For r < r where Vc OU i(^) > 2m e c 2 , we treat the term 2m e c 2 V cou \(r) in Eqs. fl6lfl"2l as a 
small correction term, and find the following inequality is always true 

< eutral (r) > nf arged (r), r < r , (15) 

where ra neutral ( r ) an d n charged^ r ^ s ^ anc [ f or electron densities of neutral and super charged 
cores. For the range r > r , nf laTSed (r) = in the case of super charged core, while 
ng 6Utral (r) — > in the case of neutral core, which should be calculated in non-relativistic 
approximation: the term V 2 oul (r) in Eq. is neglected. 
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In conclusion, the physical scenarios and Thomas-Fermi equations of neutral and super 
charged cores are slightly different. When the proton density n p of cores is about nuclear 
density, ultra-relativistic approximation applies for the Coulomb potential energy V cou \(r) ^> 
m e c 2 in < r < tq and tq > R c , and approximate equations and solutions for electron 
distributions inside and around cores are the same. As relativistic regime r ~ ro and non- 
relativistic regime r > r (only applied to neutral case) are approached, solutions in two 
cases are somewhat different, and need direct integrations. 
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FIG. 1: Potential energy-gap ±m e c 2 — V cou \(r) and electron mass-gap ±m e c 2 in the unit of m e c 2 
are plotted as a function of (r — i? c )/(10A 7r ). The potential depth inside core (r < R c ) is about pion 
mass m^c 2 S> m e c 2 and potential energy-gap and electron mass-gap are indicated. The radius r 
where electron distribution n e (ro) vanishes in super charged core case is indicated as ro— , since it 
is out of plotting range. 
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